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We describe the dynamics of three-dimensional fluid vesicles in steady shear flow in the vicinity 
of a wall. This is analyzed numerically at low Reynolds numbers using a boundary element method. 
The area-incompressible vesicle exhibits bending elasticity. Forces due to adhesion or gravity oppose 
the hydrodynamic lift force driving the vesicle away from a wall. We investigate three cases. First, 
a neutrally buoyant vesicle is placed in the vicinity of a wall which acts only as a geometrical 
constraint. We find that the lift velocity is linearly proportional to shear rate and decreases with 
increasing distance between the vesicle and the wall. Second, with a vesicle filled with a denser fluid, 
we find a stationary hovering state. We present an estimate of the viscous lift force which seems to 
agree with recent experiments of Lorz et al. [Europhys. Lett. 51, 468 (2000)]. Third, if the wall 
exerts an additional adhesive force, we investigate the dynamical unbinding transition which occurs 
at an adhesion strength linearly proportional to the shear rate. 
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I. INTRODUCTION 

We would like to understand how hydrodynamic flow 
affects vesicles and, in particular, shear flow because it 
occurs whenever a fluid flows near a surface. This is of 
relevance in modeling artificial or natural cells in mo- 
tion from the level of the dynamics of a single cell to the 
rheology of a suspension such as the blood. It will also 
enable us to understand the interplay between vesicles in 
shear flow and specific and nonspecific adhesive proper- 
ties of a membrane. This is a key feature in phenom- 
ena such as leukocyte locomotion Computational 
methods are now increasingly used to understand and 
analyze experiments studying the relationship between 
receptor-ligand functional properties and the dynamics 
of adhesion. An efficient numerical method simulating 
adhering spherical cells in flow has been used recently 
to investigate the physical kinetics of adhesion molecules 
1^. There are compelling reasons to go beyond spheri- 
cal cells, viz. (o) to understand how deformation couples 
with the dynamics of adhesion centers and (6) to study 
the interplay between the physical kinetics of adhesion 
molecules, orientation of a cell and viscous lift forces that 
act on non-spherical cells and do not act on spherical 
cells. It should however be noted that leukocyte loco- 
motion is further complicated by the role of actin poly- 
merization, ion channels and pseudopod formation which 
together play an important role in leukocyte response to 
shear stress Another practical use of this study con- 
cerns mechanosensory transduction which is the process 
by which certain cells convert fluid stresses to biochemical 
and/or electrical signals [^. The identification and ac- 



tivation of mechanoreceptors and intracellular signaling 
pathways are currently active areas of research. For an 
example of a case where our study is relevant, consider an 
experiment which investigates the influence of mechanical 
stress on the activation of specific proteins present on the 
membrane of endothelial cells. One method involves re- 
constituting the proteins into phospholipid vesicles which 
are then subjected to physiological levels of fluid shear 
stress in a viscometer Q. Quantitative predictions re- 
quire information regarding hydrodynamic stresses and 
tensions that develop in response to flow and the model 
discussed in this paper is suited for these purposes. 

We are also motivated by recent experiments on weakly 
adhering vesicles in shear flow These investigations 
reveal the lipid flow within the vesicle membrane and the 
flow field within the vesicle and near its outer surface. 
The influence of weak flow fields on the state of adhe- 
sion and the translational motion of vesicles are studied 
using bright field and reflection interference contrast mi- 
croscopy. Moreover, the hydrodynamic lift forces exerted 
on adhering vesicles are estimated and the lift forces are 
found significantly larger than estimates based on an ear- 
lier theoretical study j2j . An explanation for this discrep- 
ancy based on arguments using a lubrication analysis of 
flow between a spherical cap and a wall has been pro- 
posed recently ||]. 

Here, we show the results of a numerical simulation. 
With the above experiment in mind, we simulate an in- 
compressible vesicle in shear flow which is influenced by 
gravitational forces or nonspecific adhesion forces to be 
near a wall and by viscous lift forces to be away from the 
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wall. The main objectives of this paper are (a) to under- 
stand the very nature of viscous lift force; (6) to describe 
how the steady state shape is influenced by shear flow 
and forces due to gravity or an adhesion potential; and 
(c) to study the influence of shear in taking the system 
away from equilibrium. The numerical method simulates 
a realistic model of a vesicle. The tensions that develop 
in response to flow are determined self-consistently. 

Before describing the study, we take stock of what is 
known about the dynamics of vesicles, and about tech- 
niques that are relevant. The dynamics of vesicles belong 
to a broad class of problems called free-interface prob- 
lems which describe the dynamics of a particle consisting 
of a membrane that encloses a drop of an incompressible 
Newtonian fluid often called a capsule in the literature 

. The first problem discussed in this class was the be- 
haviour of a fluid drop in shear flow, when the effects of 
fluid viscosity and interfacial tension are taken into ac- 
count [0. In fact, much of the later work on elastic cap- 
sules or vesicles employed techniques and analysis used 
to understand the dynamics of liquid drops. The first 
two lines of attack involved analytical methods based on 
perturbation schemes for shapes close to a sphere, and 
theories using slender-body dynamics. Considerable ef- 
fort has been made in fluid drop dynamics ranging from 
drop breakup to physical influences such as surfactants 
and complex flows, for reviews see |^,^. An important 
step in understanding this nonlinear problem of three- 
dimensional deformation in external flows without mak- 
ing drastic approximations involves numerical methods 
using boundary element methods . The use of an un- 
structured grid of triangular elements to represent the 
interface has been shown to be efficient and numerically 
stable in simulating the dynamics of three-dimensional 
liquid droplets in shear flow in an unbounded fluid and 
in the proximity of a bounding plane wall This has 
been also used in simulating capsules with elastic mem- 
branes |l5| . We now compare this scene of activity with 
that regarding the dynamics of vesicles. 

In the last three decades, experimental and theoretical 
studies of the behaviour of free and bound fluid vesicles 
in equilibrium have met with great success But, the 
hydrodynamics of vesicles still poses unsolved problems. 
For example, consider a vesicle in unbounded shear flow. 
A minimal model of a vesicle involves bending elasticity 
and the constraints of fixed area and volume. Using this 
model, the dynamics of free vesicles in shear flow have 
been studied using boundary element methods Jl7| , p^ . 
Apart from features similar to liquid drops such as steady 
revolution of the surface, called tank-treading in the case 
of capsules, and a steady tilt, the numerical work sug- 
gested that shear is a singular perturbation even though 
one would naively expect that at very small shear rates 
the equilibrium shape is recovered. Numerical limitation 
prevented a conclusive examination of this regime. 

In two-dimensions, the dynamics of adhering vesicles 



has been investigated in models for chemotaxis [|9| and in 
shear flow [ pO[ using a boundary element method. Prac- 
tical considerations and a need to compare theoretical 
models with experiments call for an investigation of how 
hydrodynamic flows affect three-dimensional vesicles. 

In the next section, the physical model of a vesicle and 
the mathematical formulation of the hydrodynamics of 
vesicles including the implementation of the boundary 
element method will be briefly discussed. Then, we pro- 
ceed to three sections concerning the dynamics of vesicles 
in shear flow in the proximity of a wall. In Section III, we 
study the nature of viscous hydrodynamic lift force and 
we assume that the vesicle does not experience gravita- 
tional force, i.e., the vesicle is neutrally buoyant, and the 
vesicle is not influenced by forces due to an adhesion po- 
tential. In Section IV, gravitational forces are introduced 
and we compare our results with the experiment. In Sec- 
tion V, the influence of a nonspecific adhesion potential is 
considered. And in Section VI, our conclusions regarding 
the dynamics of vesicles near a wall are summarized. 



II. MODEL 

Consider the vesicle to be a two-dimensional surface 
embedded in three-dimensional space. The instanta- 
neous membrane configuration R(si, S2) is parametrized 
by internal coordinates (si,S2)- The energy, E, of the 
vesicle - with an area element dS - 



E = ids 



(1) 



has three contributions. The first due to the squared 
mean curvature, , describes the bending energy |l^ ] 
with bending rigidity k. The second term is due to a 
locally varying isotropic tension E which is needed to en- 
sure local incompressibility of the membrane. The third 
term is the adhesion energy due to the proximity of a 
homogeneous substrate which exerts a nonspecific inter- 
action. The geometry is schematically represented in Fig- 
ure 1. The adhesion potential, with the wall in the plane 
2; = 0, is chosen to be 



W(z) = W{dolzf [{d^/zf - 2] . 



(2) 



Here, W is the adhesion strength. The potential is re- 
pulsive as l/z*^ for z ^ do and attractive as — at 
long distance. The potential has a minimum oi — W at 
z — do. In the case of free vesicles, the adhesion term is 
absent. The membrane force density, f , reads 



f(R)^- — ^ 



1 5E 



(3) 



where g is the determinant of the metric tensor. 

To describe the hydrodynamics of vesicles, we assume 
that the Reynolds number of the flow inside and around 
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the vesicle based upon the vesicle size (say the radius of 
the spherical vesicle of equal volume) is sufficiently small 
so that the velocity, v, of the fluid and the pressure field, 
p, is governed by the equations of Stokes flow or creeping 
flow equations. Thus 



= 0, 
0, 



(4) 



where and p™ are the pressures, and 77°"* and 77™ are 
the viscosities associated with the outer and inner fluids 
respectively, and di denotes differentiation with respect 
to the coordinate Xi — x,y, z for i = 1,2,3 respectively. 
The summation convention is used over doubly occurring 
indices. 

If the vesicle is not neutrally buoyant, then we account 
for the body force acting on the vesicle due to gravity. 
This is done by modifying the membrane force density. 
The modified force density, f™°'', is given by [O] 



fmod^j^) = f (R) + _ (g . R) n(R). 



(5) 



Here, p™ and are the densities of the outer and in- 
ner fluids respectively; n is the unit vector normal to 
the vesicle, g = —go^z is the acceleration due to gravity 
where 50 — 9.81 m s~^ and is the unit vector in the 
z-direction. 

To express the no-slip condition on bounding surfaces 
and also the impermeability of the membrane, we require 
that the velocity be continuous across the vesicle. This 
also provides the kinematic condition by which the vesi- 
cle shape changes with time. The hydrodynamic surface 
force is allowed to undergo a discontinuity that is bal- 
anced by the membrane forces, f™°'*(R). 

Rather than solve for the fluid velocity at all points 
in space, it is advantageous to use a boundary-integral 
method by which the Stokes equations inside and outside 
are cast into an integral form that involves only quanti- 
ties evaluated on the vesicle surface . This formalism 
provides us with an integral equation for the membrane 
velocity in terms of f , 

^;,(R") = ur(R") + ^ ^G,;,(R",R'^)/,(R'^)d5, (6) 

where v°°{x,y,z) = (72,0,0) is the incident shear flow 
(refer Figure 1) and Gij(R",R'') is the appropriate 
Green's function for the velocity. In the case of free vesi- 
cles in unbounded space, the free-space Green's function 
is the Stokeslet; and, in the case of bound vesicles, we 
choose the semi-infinite space Green's function for the 
velocity |^l|. To make computations tractable, we as- 
sume that the viscosity of the fluid inside the vesicle is 



equal to that of the suspending fluid, (77° 



rj 



If the viscosities are different, equation (6) is modified to 
include an additional term which turns it into an inte- 
gral equation in v. Then, the equation has to be solved 



iteratively or v must be determined by matrix inversion 
at any instant of time. We also note in passing that by 
choosing the same viscosity, we do not allow for tumbling 
motion in shear flow which occurs when the dynamics is 
vorticity-dominated. For liquid drops, tumbling is shown 
to occur for A = 77™/?^°"* > 4 [|l4|. 

We characterize vesicles by a dimensionless number, 
the reduced volume 



V = y/(47ri?3/3), 



(7) 



where V is the enclosed volume and the surface area A de- 
termines the length scale i?o = A/ (47r). For a sphere, 
V — 1. The bending rigidity of the membrane k sets 
the energy scale. We do not incorporate thermal fluc- 
tuations. Throughout this paper, unless otherwise men- 
tioned, lengths are expressed in units of i?Oj the adhesion 
strength in units of k/Rq, and the shear rate in units 
of K/(87r?7i?Q). When gravitational forces are considered, 
the dimensionless gravity parameter 



g' ^ ir - p'"'')goRt/^ 



(8) 



measures this effect. In simulations involving the adhe- 
sion potential equation (2), we take do — 0.01i?o- 

The implementation of the boundary element method 
is now well-established p^ , p^ 171^]. Here, we briefly 
mention the essential steps. 

(a) The surface is approximated by a grid of triangular 
elements. We take 512 triangles with 258 nodes. 

(b) The force density is calculated at the nodes. 

(c) Along with the condition dictated by 2D incom- 
pressibility. 



d_ 

dt 



y/g = diVi - nitijdjV^ = 0, 



(9) 



the integral equation equation (6) yield equations for the 
local tensions S. Unlike droplets with constant tension, 
the tensions S develop so that the membrane element 
deforms while maintaining the original area. 

(d) After computing S, the integral equation (6) is 
solved for v. 

(e) The position of the nodes are advanced, the coordi- 
nates updated and we return to step (b) unless a steady 
state is reached. 

Before tackling cases involving externally imposed hy- 
drodynamic flows, we test if the numerical scheme yields 
the stationary shapes that are known to exist in equi- 
librium conditions, i.e., neutrally buoyant free vesicles 
JlGf . The minimal model we have chosen allows a wide 
spectrum of shapes. For 0.65 ^ v ^ 1, the stable states 
are prolate-like spheroids. For v < 0.75, oblate disco- 
cytes are locally stable. So far, there is no indication of 
the presence of locally stable non-axisymmetric shapes 
in this model. For v°° = 0, our numerical scheme gives 
the stable and metastable axisymmetric shapes, includ- 
ing dumb-bell shaped prolates and erythrocyte-like bi- 
concave discocytes. In Figure 2, we show an example 
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in which the initial oblate spheroid relaxes to the final 
equilibrium (stationary) biconcave shape. Within this 
minimal model, the vesicle relaxes to the next dynami- 
cally accessible locally stable shape. Depending on the 
initial condition, this is either a prolate, an oblate or a 
stomatocytc. 

In contrast to vesicles, the shape of incompressible in- 
terfaces that are not influenced by curvature elasticity is 
infinitely degenerate. In such cases, the stationary shape 
depends only upon the membrane surface area, and as 
a consequence of neglecting elastic behaviour, the inter- 
face can retain any shape with the same surface area. At 
the other end of the spectrum, liquid drops with inter- 
faces influenced by a constant surface tension and also 
compressible interfaces which are only influenced by cur- 
vature elasticity remain spherical. 



III. UNBOUND NEUTRALLY BUOYANT 
VESICLES NEAR A WALL 

We now simulate the influence of steady shear flow on 
vesicles in the proximity of a wall. To retrieve the bare 
contribution of hydrodynamics to the lift force, adhesion 
and gravity are not considered in this section. The vesi- 
cle is initially at a small distance from the wall and then, 
shear flow is imposed. In Figure 3, a typical sequence of 
"snapshots" are shown. In the far-field, the vesicle trans- 
lates with a speed of the same order of the difference of 
velocity, due to external shear flow, that exists across 
the vesicle. The vesicle also experiences a lift velocity. 
From Figure 3, it can be extracted that the lift velocity 
is an order of magnitude smaller than the translational 
velocity. 

Analyzing several cases, we find the main characteris- 
tics in this setup to be: 

(«) the membrane of the vesicle tank-treads; 

(m) the vesicle develops a steady tilt which is roughly 
independent of the shear rate; and 

{Hi) the shape changes to a prolate-like ellipsoid; 
This is similar to the dynamics of a free vesicle in un- 
bounded shear flow fl^ . The first two features are un- 
derstood The tank-treading motion is due to the 
rotational component of linear shear flow. The balance of 
the moments on the vesicle due to the effect of shear-flow 
acting on a stationary inclined vesicle and that due to 
the tank-treading motion results in a steady tilt. These 
two moments are linearly proportional to shear rate and 
hence, the tilt is independent of the shear rate in the 
steady state. In the case of liquid droplets, the tilt re- 
duces more drastically with shear rate because of the 
elongation of the droplet, and thereby a reduction in v. 
But in the case of vesicles considered here, t; is a constant 
due to incompressibility. As mentioned in the introduc- 
tion, for the third feature regarding the prolate ellipsoidal 



shape, it is still not understood whether it will persist for 
arbitrarily small shear rates. 

From the above discussion, it is clear that the steady 
tilt of the vesicle plays the leading role in breaking the 
fore-aft symmetry of the vesicle with respect to shear 
flow. Since the tilt is independent of shear rate, the ex- 
cess pressure in the space between the vesicle and the wall 
can be expected to cause a lift force which is proportional 
to shear rate, 7. In fact, the dynamics of a vesicle near a 
non-adhering wall is roughly similar for any shear rate; 
and as a function of r = jt, the "snapshots" as shown in 
Figure 3 will be similar for vesicles with the same reduced 
volume. In Figure 4, we plot for different shear rates the 
ratio of lift velocity, vu ft , to shear rate against the mean 
height. The data collapse indicates that vuft oc 7. 

We now derive an empirical expression for the lift ve- 
locity guided by the far-field result obtained for ellip- 
soidal cells: VHft = UjiR^/h^) where J7 ~ 0.1 - 0.3 for 
V ~ 0.9 — 0.99 [Q. Assuming the same power law, we 
get an estimated asymptotic fit for the lift velocity to be 

vuft^ 0.08 j-^ {v = 0.95) (10) 

Here, h is the mean distance of the center of the vesicle to 
the substrate. We have simply calculated the mean cen- 
ter and the mean velocity to be the arithmetical mean of 
the coordinates and velocity of the nodes respectively. 
Then, vuft is just the vertical or z-component of the 
mean velocity. The fit needs to be taken with caution 
due to three reasons. First, the exact dependence of un/t 
on h is not easy to access with the limited range as shown 
in Figure 4. Second, the asymptotic fit (and also, far- 
field results based on analytical methods) for vesicle mi- 
gration assumes that the problem can be treated in a 
quasistatic way. This implies that it is strictly valid only 
when the time-scale for deformation and shape changes 
is much smaller than the time-scale for migration. Since 
our numerical method is similar to an experimental situ- 
ation and we do not arbitrarily fix the position of a vesi- 
cle, the fit is probably valid only for large shear rates. 
This is partially evident from Figure 4 where we see that 
the plots for lower shear rates do not exactly collapse 
onto the plots for higher shear rates. Third, there is 
no reason to expect asymptotic results to hold good for 
h ~ 0.25 — 0.5i?o- But, there is also no reason to be 
alarmed if it is true in that range. In fact, analogous 
results for liquid drops do indeed match remarkably well 
for A = 1 at /i ~ 0.25 - 0.5i?o We can also infer 

from far-field results that the sign of the prefactor is de- 
termined by the sign of the tilt 9 of the vesicle. Hence, 
when 7 —7 and 9 —9, vuft vuft which agrees 
with expectations that the direction of migration should 
be independent of the direction of shear flow. 

We emphasize that the dominant effect is due to the 
tilt of the vesicle and this leads to a lift force, Fuft oc 7. 
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Additional effects due to deformation can be expected to 
be of 0(72). 

For h ^ Rq, the asymptotic lift velocity can be taken 
to be the lower bound for the lift velocity. The actual 
velocity could be expected to be higher due to transient 
effects or due the effect of "touching" the wall. Assum- 
ing Stokes' result Q^n-qR^vuft, this lower bound for the 
viscous lift force can be estimated from equation (10) for 
V = 0.95 to be 

Fuft ^ 0.45^777-^ (11) 

With this estimate of the lift force, a comparison can be 
made with the results of the experiment. We shall do so 
only at the end of the next section, after investigating 
the influence of gravity. 

IV. INFLUENCE OF GRAVITY 

The effect of gravity on vesicles arises from the fre- 
quently employed experimental technique to stabilize the 
vesicle at the bottom of the measurement chamber by a 
difference in density between the fluids inside and outside 
the vesicle. In the experiment [Q , this is done by choosing 
iso-osmolar but different buffers for swelling the vesicle 
and suspending the vesicle. The threshold velocity at 
which unbinding of the vesicle from the substrate occurs 
allows to estimate the lift force. After unbinding, the 
vesicle hovers at a distance at which electrostatic forces 
or van der Waals forces that cause adhesion are negligi- 
ble. The lift force in the hovering state can therefore be 
assumed to be counteracted by gravitational forces only. 

We simulate the influence of gravitational force on a 
vesicle. As in the previous section, we place the vesi- 
cle initially at a distance ~ 0.1i?o — 0.5-Ro ^nd then, 
shear flow is applied. It would of course be more realis- 
tic if we started with a vesicle which rested on the wall, 
but at gravitational strengths which are of relevance, this 
initial conflguration leads to numerical instabilities. In- 
stead, the initial is allowed to be non-zero but we check 
if the steady state dynamics is independent of /iq and the 
simulations conflrm this view. We note in passing that 
the steady state depends on the reduced volume of the 
vesicle but does not depend on the initial shape. 

As expected, gravity counteracts viscous lift force and 
the vesicle hovers at a distance away from the wall. We 
now study how this distance depends on the shear rate. 
It is clear from Figure 5(a) that the initial /iq can be cho- 
sen so that we can optimize the computational time in 
evolving the vesicle to the steady state. At large shear 
rates, the vesicle is sufficiently far from the wall and as 
discussed in the previous section, the lift is determined 
by the tilt alone. Meanwhile, at low shear rates, the vesi- 
cle is also deformed by the proximity to the wall. The 
dependence of the steady state mean height on shear rate 



is shown in Figure 5(b) and seems to indicate those two 
regimes. 

Let us now consider the three-dimensional shape of the 
vesicle in the steady state. In Figure 6, we show a com- 
parison of the profile and also the top view of the vesicle 
for different shear rates. We note the features that are 
immediately apparent from Figure 6. 

[i) Vesicles with small excess areas w 1 do not ex- 
hibit large deformations with increasing shear rate. 

[ii) At large shear rates [refer Figures 6(c) and 6(d)], 
the shape of the vesicle is roughly independent of shear 
rate. 

[Hi) At any shear rate, the part of the vesicle away 
from the wall shows roughly the same inclination. 

{iv) At lower shear rates [refer Figures 6(a) and 6(b)], 
gravity reduces the vesicle-wall distance. To allow for 
greater contact area, the shape changes from a prolate- 
like to an oblate-like ellipsoid. 

{v) The tilt of the contact area decreases with decreas- 
ing shear rate. There is no unique definition of this tilt 
and therefore, the tilt of the "base" is difficult to quantify. 

Now, we compare our results with the experiment [Q. 
There are some problems in making this comparison. The 
reduced volume v and the gravitational strength g' are 
not exactly known. For vesicles with very small excess 
area, the experiment does not report any measurable tilt 
of the contact area nor the overall shape. Such a picture 
fits in with our description only in the limit — > 00. 
We take the shape of the vesicle to be close to a sphere 
with V = 0.95. We assume the values quoted, viz. the 
density difference, p™ — ~ 5.2 mg/ml, the radius of 
the vesicle i?o — 14^m, the vesicle-wall gap ~ lOOnm 
and shear stresses typically 777 ~ 1.2mPa. These experi- 
mental values lead to an estimated gravitational force on 
a vesicle with volume V of 

Fgrav = (p™ - P°"*)5o^ - 6 X lO-i^N. (12) 

As a theoretical estimate, we obtain the lift force based 
on equation (11) assuming that /i ~ i?o as 

F,,/t~ 0.457r777i?2~3xl0-i3N. (13) 

The agreement should be taken with a grain of salt if one 
takes into consideration all the assumptions that have 
gone into these estimates. 

V. INFLUENCE OF ADHESION 

Now, we consider the dynamics of vesicles adhering to 
a surface due to nonspecific interactions. Here, the inter- 
play between forces due to an adhesion potential, equa- 
tion (2), and hydrodynamic forces is considered. The 
vesicle is assumed to be neutrally buoyant. 

Like an experiment, we start with an adhering vesicle, 
then apply shear flow and watch how the vesicle responds. 
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The dynamical evolution in the transient stage is shown 
in Figure 7. These "snapshots" are qualitatively similar 
for any adhering vesicle with arbitrary reduced volume 
at any adhesion strength or shear rate. The tilting of the 
vesicle is similar to the scenario in the previous section. 

The steady state of the adhering vesicle, of course, 
depends on the reduced volume, adhesion strength and 
shear rate. Unlike the case involving gravity, vesicles ad- 
hering to a substrate by the influence of a short-range po- 
tential cannot be made to hover at any arbitrary height 
by tuning the strength of adhesion. There is a critical 
shear rate above which the vesicle unbinds from the wall. 
For smaller shear rates, the vesicle tank-treads along the 
substrate. In Figure 8, we show how the critical shear 
rate varies with increasing adhesion strength for a vesicle 
with reduced volume, v = 0.95. For 7 ^ 7c, the steady 
state shape of the vesicle resembles the shape shown in 
Figure 7(c). Two main features of this state are to be 
noted. One, the tilt of the vesicle is nearly the same as 
that of a non-adhering vesicle near a wall, and not very 
different from the tilt of a free vesicle in shear flow. Two, 
the surface of adhesion is nearly zero. The vesicle can 
be said to be "pinned" . This state is nearly similar for 
any adhesion strength and this explains the roughly lin- 
ear relationship between the critical shear rate, 7c, and 
the adhesion strength. For lower shear rates or higher 
adhesion strengths, the area of contact increases and a 
typical example of a "bound" steady state is shown in 
Figure 9. We find that the process of unbinding, viz. 
"bound" — > "pinned" — *■ unbound or "free" , is similar to 
that described with simulations of the two-dimensional 
model @. 

Finally, we summarize the results of the simulations 
with a dynamical phase diagram. First, we recall the 
behaviour of adhering vesicles under equilibrium condi- 
tions. In the absence of shear, it is known that there are 
two states jl^. For weak adhesion, the shape of the vesi- 
cle resembles the free shape and the surface of adhesion 
is zero. This state is called the pinned state. It should 
however be noted that there is no unique definition of 
the pinned state in the dynamical picture. When adhe- 
sion is strong, there is the bound state with non-zero 
adhesion area. The adhesion area also depends on the 
excess area which is a measure of the reduced volume. 
For 0.5 < f < 1, the transition from the bound axisym- 
metric shape to the pinned state is discontinuous. Keep- 
ing in mind the similarity between the equilibrium and 
dynamical situations, we can construct a schematic dy- 
namical phase diagram (refer Figure 10). Apart from the 
two-stage dynamical unbinding, we expect in the limit 
of weak adhesion, a regime in which the unbinding goes 
through only one stage. In the case of an adhesion poten- 
tial with finite range, there will be a finite critical shear 
rate for unbinding a vesicle which is pinned under equi- 
librium conditions. If the vesicle is pinned by a contact 
potential, the vesicle will unbind for any shear rate. This 



regime is difficult to realize either numerically or experi- 
mentally. In principle, there could be a direct transition 
between bound and free vesicles but we have not found 
any evidence for it based on the numerical simulations 
over a wide range of parameters. 

VI. CONCLUSION 

We conclude this study of the dynamics of vesicles near 
a wall with observations regarding the computation, sum- 
mary of main results, and certain remarks about current 
limitations and ongoing research. Using the boundary el- 
ement method with a grid of triangular elements, we are 
able to simulate incompressible vesicles in bounded and 
unbounded shear flow. With moderate shear rates, ad- 
hesion and gravitation strengths, the numerical method 
is stable and effective in studying the steady state of ad- 
hering, hovering or free vesicles. The algorithm which is 
used to ensure area incompressibility is the main hurdle 
in improving the speed of computation and also in using 
finer discretization. 

We are able to find the following features regarding 
the nature of viscous lift force on incompressible vesicles 
with bending elasticity: 

{i) The lift force due to the tilt of the vesicle is linearly 
proportional to the shear rate. It is the dominant part 
in our simulations also for small shear rates. 

(m) Our numerical results agree well with a recent ex- 
periment. 

[Hi) The general features of dynamical unbinding is 
similar to that in simulations of the two-dimensional 
model even though the additional possibility of a prolate- 
oblate shape change complicates the issue. 

We also remark about interesting issues that need to 
be resolved, both experimentally and theoretically. We 
have only considered a situation where the viscosities of 
the inside and outside fluid are equal. This is definitely 
an oversimplification of experiments with cells though 
it is not a limitation for experiments with phospholipid 
vesicles. We still have to understand whether the model 
considered is singular in the limit of very small shear 
rates. We have not done simulations with the ratio of 
vesicle-wall gap to radius, ft,o/-Ro < 0-001. In this limit 
of vanishing gap, features such as the roughness of the 
wall and also, the no-slip boundary condition will have 
to be reassessed. Finally, we note that the numerical 
method discussed here is amenable to an extension that 
includes specific receptor-ligand type of adhesion. 
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FIG. 1. Schematic representation of the geometry show- 
ing a vesicle, the wall in the plane z = and applied linear 
shear flow. On extreme right, the mean center {xc,0,zc) is 
shown as C. The height of the center, h, and the vesicle-wall 
distance, ho, are also shown. 
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(a) 



FIG. 2. An example of a test of the numerical method in 
which an initial shape relaxes to a stationary shape known to 
exist in the absence of externally applied flow, (a) The initial 
shape is chosen to be an oblate convex spheroid with reduced 
volume, V = 0.6. (b) After relaxation, the equilibrium or 
stationary shape is a stable oblate biconcave discocyte with 
the initial area and volume. 
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FIG. 3. "Snapshots" of the hft-ofF of an unbound vesi- 
cle away from a waU in the presence of shear flow. The 
plane z = shown is the wall. The x-coordinate at the 
endpoints of the wall is shown at each instance of time, t. 
Here, v = 0.95, 7 = 30. (a) At t = 0, the initial shape is an 
oblate-like spheroid. The vesicle is at a distance ho = 0.1-Ro. 
(b) The shape is now a prolate-like ellipsoid and the vesicle 
is tilted with respect to the shear plane. Here, t = 0.075. (c) 
The shape and tilt of the vesicle is roughly the same as before, 
and the lift-off is now clearly seen at t = 1.5. All quantities 
are given in dimensionless units as discussed in the text. 
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FIG. 4. The ratio of mean lift velocity, vaft, to shear rate 
times the radius of vesicle, 7-Ro, is plotted against the ratio of 
the mean height h of the vesicle from the wall to the the radius 
of vesicle for different shear rates (definitions are given in the 
text). The lift velocity of the vesicle is roughly proportional 
to the shear rate. Here, v — 0.95, and in dimensionless units, 
7 — 6(o), 10(+), 30(Q), 50(x) and the estimated asymptotic 
fit given by equation (10) is plotted using the asterisk (*). 
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FIG. 5. (a) Mean vertical height of the vesicle h versus 
for different shear rates. The shear rates, 7, are (a) 64 (6) 96 
(c) 150 (d) 200 (e) 300 (/) 350 (g) 400 (h) 500. (b) Steady 
state height of the vesicle h versus the shear rate, 7. Here, 
V — 0.95 and 7 is given in dimensionless units as described in 
the text. 
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FIG. 6. The profile of a vesicle in the x — z plane and the 
top view (looking down z-axis) are shown for four diflterent 
shear rates (in diinensionless units), (a) 7 = 64 (o), (b) 7 = 96 
(+), (c) 7 = 350 (x) and (d) 7 = 400 (□). The profile is 
shown on axes Xc — Zc such that their mean centers coincide. 
Here, v = 0.95 and g — 6.4. 
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X = —1.5 X = 1.5 




FIG. 7. The dynamical evolution of an adhering vesicle is 
shown. Here, v = 0.95, 7 = 2 and W = 0.1 (all quantities 
scaled to be dimensionless as described in the text). The 
top and side views are shown on the right side and left side 
respectively, (a) An adhering vesicle in equilibrium resembles 
an oblate spheroid at t = 0. (b) In the transient stage the 
shape changes to a tilted prolate-like ellipsoid, [t = 0.6). 
(c) In the steady state, t = 1.2, the tank-treading vesicle 
is "pinned" to the wall and slips/rolls along the wall. The 
wall in the plane z = is in the side view. The a;-coordinate 
at the endpoints of the wall is shown at each instance of time. 
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FIG. 8. The critical shear for unbinding, 7c, is plot- 
ted against the adhesion strength W for v = 0.95, and 7c 
is linearly proportional to the adhesion strength W . Here, 
do = 0.006. With K = lO'^V, r) = 10"^ Js/m^ i?o = S^tm, 
and with W roug hly 10^"' - 3 X 10"'* J/m^, the critical shear 
rate is roughly in the range 0.1 — 4s ~^ 
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FIG. 10. Schematic phase diagram in the shear rate 7, ad- 
hesion strength W plane with free, pinned and bound shapes 
for vesicles with arbitrary area and volume. 
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